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Abstract 

We describe a numerical procedure to resum multiple emission effects in event 
shape variables and jet rates. 

1 Introduction 

Event shape variables (and jet rates) are among the most studied observables 
in QCD. They are useful both to measure as 0] and to search for genuine non 
perturbative (NP) effects Q. A perturbative (PT) fixed order expansion describes 
well an event shape rate (the fraction '^{v) of events whose value of the event shape 
V is less than v) in the region t; ~ 1. In the less inclusive region f -C 1 large 
logarithms a™ In" v arising from incomplete real-virtual cancellations have to be 
resummed to all orders to give meaning to the PT series [3-7]. In order to fix the 
scale of Os one needs to control in In S all terms of the form a" L""*"^ (leading LL or 
double DL logarithms) and a^L" (next-to-leading NLL or single SL logarithms), 
with L = \nl/v. It is therefore important to understand the sources of large 
logarithmic contributions and to find a general way to resum them, at least at 
NLL level. The aim of the paper is to describe a general procedure to resum SL 
terms arising from multiple emission effects, which is the most difficult task in 
any resummation programme. 



2 Classification of large logarithms 

Leading logarithms originate from soft and collinear gluon radiation and are in 
general straightforward to resum. SL contributions instead have a variety of 
sources: 

• running of the coupling; 



• soft emission at large angles; 

• hard collinear splitting; 

• multiple emission effects. 

The first three contributions can be easily resummed by computing the first or- 
der result (in the soft-collinear limit), and exponentiating the answer. The last 
contribution is the most difficult to address. In the luckiest case (additive observ- 
able) it requires only one Laplace transform in more complicated cases one 
needs to introduce an additional amount of Fourier transforms (which can be 
as many as five in the thrust minor case |^]). There are situations (thrust major, 
oblateness) in which an analytical formulation does not even exist at all. It is 
therefore important to investigate if there is a general method capable to resum 
SL terms arising from multiple emission effects. 



3 Resummation of multiple emission effects 

A resummed integrated distribution for a 'suitable' (for a definition see §) ob- 
servable V can be written in the form 

j:(v) = e-«(^) HR') , R' = -v^ , (1) 

dv 

where the radiator R builds up the Sudakov form factor obtained by exponentiat- 
ing the one soft-collinear gluon contribution, while J~{R') is a SL function which 
accounts for multiple emission effects. Suppose that we know the resummed 
distribution Ds{vs) of a 'simple' variable Vg- The distribution D{v) of a more 
complicated observable V is related to Ds{vs) by: 

D{v) = V [ ^P{v\vs) DM , (2) 

with P(v\vs) the probability of having V equal to v given a value Vg for Vg. The 
result (|2|) is quite general. If V and Vg have the same DL structure, it can be 
simplified to give, at SL accuracy: 

^(v) = ^s{v) T{R') , T{R') = I ^x-%(x, R') , p (f , = ^P{v\vs) ■ 

(3) 

If we now choose Vs{{ki}) = maxj{y(fcj)}, the resummation of Ss(^^s) leads di- 
rectly to 

^s{vs) = e-^^''^) , (4) 

so that the function J-{R') obtained from (^) coincides with the one introduced 
in (0). 

The effort of any resummation is then to compute p(x, R'). This can be done 
numerically with the following Monte Carlo procedure. One starts with a given 
Born configuration with a fixed number of jets. One then 

fixes a value Vs of the simple observable; 

1 generates a soft and collinear emission according to the phase space R' , with 
V{ki) < Vih.i); 



II if V{ki) < e ^ 1 one stops, otherwise goes back to step /; 

III given the momentum set {A;j}j=i_...^„ one computes the value of the observ- 
able V = V{ki, . . . , kn)- 

The above procedure gives the probability P{v\vs), which can be integrated ac- 
cording to @ to give the function J^{R'). 
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Figure 1: Comparison of the analytical and numerical results for the function J-'{R') 
for the thrust and the two broadenings. The lines represent the analytical predictions 
while the dots are the numerical results. 

The method is process independent and has been tested in the case of observ- 
ables for which an analytical resummation is known. For the thrust and the two 
broadenings Bt and B\y the results are shown in figure ||. The Monte Carlo pro- 
cedure reproduces with great precision the function J-{R') for these observables. 
This means that numerical results are truly interchangeable with analytical ones. 
The method has then been applied |8| to the thrust major Tm and the oblateness 
O, for which a resummed prediction did not exist yet, and to the three-jet reso- 
lution ^3 (Durham algorithm pl), for which only a part of SL contributions had 
been computed ||To[| (see figure^2[). This new method can then be seen as a first 
step towards a full automatic resummation of event shapes and jet rates at NLL 
accuracy. 
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Figure 2: The function jF(i?') for the thrust major Tm, the oblateness O and the 
three-jet resolution y^. 
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